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Abstract 

We examine the prospects for determining tan (5 from heavy Higgs scalar pro- 
duction in the minimal super symmetric standard model at a future e + e~ collider. 
Our analysis is independent of assumptions of parameter unification, and we con- 
sider general radiative corrections in the Higgs sector. Bounds are presented for 
yfs = 500 GeV and 1 TeV, several Higgs masses, and a variety of integrated lumi- 
nosities. For all cases considered, it is possible to distinguish low, moderate, and 
high tan (3. In addition, we find stringent constraints for 3 ^5 tan ^ 10, and, for 
some scenarios, also interesting bounds on high tan/3 through tbH^ production. 
Such measurements may provide strong tests of the Yukawa unifications in grand 
unified theories and make possible highly precise determinations of soft SUSY 
breaking mass parameters. 
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1 INTRODUCTION 



Supersymmetry (SUSY) is an attractive target of future high energy experiments, and 
the discovery of supersymmetric particles is eagerly anticipated at proposed new collid- 
ers. The discovery of superparticles is, however, not the end of the story. Once SUSY is 
found, we will be at a new stage in high energy physics. First of all, we should check the 
supersymmetric relations among various parameters to really confirm supersymmetry. 
We can also begin the exciting investigation of physics of high energy scales by using 
renormalization group analysis: the SUSY breaking parameters as well as the coupling 
constants contain information about high energy scales, and they can give us some hints 
about new physics, such as grand unification (GUT), flavor symmetry, or the mechanism 
of SUSY breaking. For these programs, accurate determinations of the parameters in 
the lagrangian are crucial and essential. 

Among the various parameters in SUSY models, tan/3 is one of the most important. 
One reason is that tan/3 plays a significant role in relating experimental observables 
to the parameters in the lagrangian. For example, reconstructions of the Yukawa cou- 
pling constants, mass matrices of the charginos and neutralinos, and the SUSY breaking 
masses of sfermions require a knowledge of tan (3. Furthermore, a precise determination 
of tan (3 can be a check of various models that prefer a specific value of tan (3. 

In this study, we discuss the prospects for tan (3 determination from the production 
and decay of Higgs bosons at the Next e + e~ Linear Collider (NLC) ^ 0. Detailed study 
of Higgs properties can give us a good determination of tan/3, since the interactions of 
the (heavy) Higgs bosons with quarks and leptons depend on tan/3. This possibility 
has been studied in a model independent approach || and in the framework of the 
minimal supergravity J4|. The NLC will be a good place for the detailed study of Higgs 
properties, and ~ 1000 heavy Higgs pairs can be produced if the heavy Higgs bosons 
are kinematically accessible. By using the standard collider parameters (a/s = 500 GeV 
and L = 50 fb _1 /yr for phase I, and y/s = 1 TeV and C = 200 fb _1 /yr for phase II), we 
estimate the expected accuracy of the measured tan (3 as a function of the actual value 
of tan/3. We will see that the error can be as small as O(10 %) or less for most of the 
theoretically interesting tan (3 values. 

In determining tan/3 with Higgs bosons, the model dependence is weak, i.e., to deter- 
mine tan (3 from Higgs bosons, we require only that (heavy) Higgs bosons be produced 
at the NLC. In order to make this point clear, we do not make any assumption that 
strongly depends on some specific model, such as the minimal supergravity model. In 
our analysis, we assume only that nature is described by the field content of the minimal 
supersymmetric standard model (MSSM). All the relevant parameters which are needed 
to determine tan/3 are then required to be measured experimentally. Furthermore, we 
emphasize that our approach results in quite an accurate measurement of tan/3 if the 
actual tan/3 is in the moderate region (~ 3 - 10). Several other methods have been 



1 



proposed to determine tan/3, such as those using charginos 0, staus ||, or the muon 
(g — 2) |7|. The discovery of H, A — > rf at the LHC may also be used to set a lower 
bound on tan (3 |§ , though, in general, the study of the heavy Higgs sector appears to 
be one of the most challenging for the LHC 0. However, all of these methods do not 
give good results if the underlying tan (3 is in the moderate region. Thus, our method is 
complementary to the other analyses. 



2 HIGGS BOSONS IN THE MSSM 



First, let us briefly review the Higgs sector in the MSSM ||10|| . The MSSM contains two 
Higgs doublets: 

When the neutral components of these Higgs fields obtain vacuum expectation values 
(VEVs), electroweak symmetry is broken. One combination of VEVs is constrained 
so that we obtain the correct value of the Fermi constant: 2((if°) 2 + (H®) 2 ) = v 2 ~ 
(246 GeV) 2 . On the other hand, their ratio is the free parameter which is tan/?: 

tan/3 =(H° 2 )/(H° 1 ). (2) 

By expanding the Higgs fields around their VEVs, we obtain physical Higgses as 
well as the Nambu-Goldstone bosons. In order to obtain the physical modes, it is more 
convenient to use another basis $i and $2: 





(3) 



cos (3 sin (3 
sin/3 cos/5 j y iu ±± 2 

In this basis, $1 gets a VEV, while $2 does not. We expand $1 and $2 as 

ti = /(, + A-HG°)/V2\ (4) 

* 2 = ( ^- + ' H A y^ ) . (5) 

Then, from the fact that electroweak symmetry is broken by the VEV of $1, the Nambu- 
Goldstone bosons G° and G ± are contained only in $ x . The other fields <p2, A, and 
are physical degrees of freedom. The pseudoscalar A and charged Higgs H ± are 
mass eigenstates. On the other hand, the CP-even scalars, </>i and 02, mix in the mass 
matrix. Mass eigenstates, h and H, can be obtained by a unitary transformation: 




sin(/3 — a) cos(/3 — a 
- cos(/5 — a) sin(/3 — a 




(6) 
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where the unitary matrix is parametrized by a new parameter a. We define h to be 
lighter than H. 

At tree level, mass and mixing parameters are related to each other; once we fix ttla 
and tan /3, all the masses and the mixing parameter a are fixed. However, the tree level 



relations can be significantly modified by radiative corrections [II], and hence it may 



be dangerous to assume tree level relations in the analysis. Therefore, we regard all the 
masses and mixings as parameters to be measured by experiments, and uncertainties in 
these measurements enter our analysis as systematic errors. 

Here, we comment on the so-called "decoupling limit" of the heavy Higgses. When 
rriA is much larger than rrtz, the mixing between $1 and $2 becomes small: cos((3 — a) — > 
0. In this limit, $1 behaves like the standard model Higgs, while the heavy Higgses (H, 
A, and H^ is like an extra doublet with degenerate mass. In our study, we assume 
that the charged Higgs mass is heavier than the top quark mass, which implies that the 
decoupling limit is more or less realized. 

In the decoupling limit, h is mainly produced in association with the Z boson 
(e + e _ — > Zh), while for the heavy Higgses, pair productions (e + e~ — > AH, e + e~ — > 
H + H~) are the most important processes. The cross section for H + H~ production is 
independent of a and /3, while those for Zh and AH are both proportional to sin 2 (/3 — a). 
From the precise measurement of the cross section of the process e + e~ — > Zh, sm 2 (fl — a) 
is well determined with accuracy ~ 2 % 0. Therefore, the cross sections can be esti- 
mated with small errors in this study. Notice that the cross sections for other processes 
(e + e~ — > ZH, e + e~ — > Ah) are also calculable, but they are too suppressed to be 
important since they are proportional to cos 2 (/3 — a). 

The Higgs bosons Hi and H 2 are responsible also for fermion masses. They have 
interactions of the form 



L Y = y t H 2 q L t c R + y h H x q L b c R + y T H x l L r c R , (7) 

where the y t , yb, and y T terms are the Yukawa couplings for m t , m b , and m T , respectively^. 
By substituting the mass eigenstates into H% and H 2 , we obtain the interactions of the 
physical Higgs bosons with the fermions. As noted before, h behaves like the standard 
model Higgs in the decoupling limit, so its interactions are insensitive to tan/3. On 
the other hand, the interactions of the heavy Higgses (H, A, and H^) with fermions 
strongly depend on tan/3. For example, the coupling of H ± to t c R and is proportional 
to cot (3, while that to t L and b c R and that to r R and v T are proportional to tan/3. Thus, 
if tan (3 is small, charged Higgs bosons mainly decay into top and bottom quarks, but 
the branching ratio of H decaying into r and v T increases as tan /3 gets large. Similarly, 

*The interactions of the Higgs bosons with the first and second generations are too suppressed to be 
important in our analysis, and we neglect them. 
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branching ratios of H and A are also sensitive to tan/3|j]. Therefore, if we measure the 
branching ratio of the heavy Higgs bosons, we can constrain tan (3. 



3 EXPERIMENTAL SIMULATION 

In this section, we present our basic idea for determining tan (3. As stated in the previous 
section, the interactions of the heavy Higgs bosons depend on tan/3, and hence the 
branching ratios of the Higgses are sensitive to tan (3. This implies that measurements 
of cross sections for heavy Higgs production with various final states can give us some 
information about tan j3. 

For this purpose, we use the following observations. First of all, Higgs bosons mostly 
decay into particles in the third generation. Thus, a large number of 6-jets is expected 
in Higgs production events if the Higgses decay into hadrons. In the NLC, 6-jets are 
expected to be selected with a high 6-tagging efficiency (e^ ~ 60 %) [0, so this can 
effectively reduce the background. Furthermore, if the charged Higgs decays into rz/ T , a 
single high energy lepton (or energetic hadrons with very low multiplicity) is expected as 
a decay product of r. This can be a striking signal of charged Higgs production followed 
by leptonic decay. The last point concerns the large tan/3 limit. If tan/3 becomes larger 
than ~ 10, all the branching ratios of the heavy Higgses lose their sensitivities to tan /3, 
and tan (3 cannot be constrained well from pair production processes. However, in this 
case, the cross section of the process e + e~ — > tbH^ can be enhanced enough to be 
observed, since the tbH^ vertex is proportional to tan/3 in the large tan/3 limit. Thus, 
in this region, the t6if ± production process becomes useful. 

Based on the above arguments, we use the following types of channels in our analysis: 

1 : 2b + / + g's+kinematical cuts to select the U H + H~" mode. 

2 : 2b + / + g's+kinematical cuts to select the "tbH ±n mode. 

3 : 36+1/ (+g's). 

4 : 36 + 0,2,3,... I (+g's). 

5 : 46. 

6 : 46+1/ (+g's). 

7 : 46 + 0, 2, 3, ... / (+g's) (but not 46). 
8: 56(+/ + g's). 



tin fact, the Yukawa couplings of H depend on a. In our numerical study, we include this correction. 
The a dependence of the interaction becomes weak in the decoupling limit. 
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In this list, "6" and "g" denote hadronic jets with and without a 6-tag, respectively, "P 
denotes an isolated, energetic e, /x, or r, and particles enclosed in parentheses are op- 
tional. In our analysis, we assume that hadronically-decaying r leptons may be identified 
as leptons, ignoring the slight degradation in statistics from multi-prong r decays. 

Channel 1 is intended to select charged Higgs pair production events with tbrv T final 
states, while channel 2 is designed to isolate ibH^ 1 production with H ± — > ru T . These 
two channels have the same event topology. Furthermore, top quark pair production 
may contribute to these channels as a significant background, if followed by the decays 
ti — > (bW + )(bW~) — > (blvi)(bqq') . It is crucial to distinguish these processes, so we 
impose certain sets of kinematical cuts for this purpose. In imposing the kinematical 
cuts, the important point is that, in these events (channels 1, 2 and the background), 
one top quark decays only into hadrons, so we can reconstruct the top quark system 
from the hadrons after first judiciously choosing one of the two 6-jets [H]. In the ti pair 
production case, the energy of the top quark is, in principle, equal to the beam energy, 
while for channels 1 and 2, it tends to be smaller than the beam energy. Thus, we 
demand the reconstructed top energy to be well below the beam energy to eliminate the 
ti background. In order to distinguish channels 1 and 2, we check the total energy of 
the hadrons. In channel 1, all the hadrons are the decay products of one charged Higgs, 
and hence the the total energy of the hadrons is equal to the beam energy. On the other 
hand, in channel 2, the total energy of the hadrons is likely to be larger than the beam 
energy. Based on these two observations, we impose kinematical cuts on channels 1 
and 2. Such kinematical cuts are highly effective once the relevant cut parameters are 
optimized, and channels 1 and 2 can be well separated. (For details, see Ref. ||.) 

Channels 3-8 receive contributions mainly from Higgs pair production events fol- 
lowed by hadronic decays of the Higgses. In these channels, we choose only events with 
large numbers of 6-jets. The resulting backgrounds have been calculated in Ref. and 
are found to be sufficiently suppressed. Therefore, we do not impose any kinematical 
cuts on these channels. 

Once the relevant channels are chosen, we can quantitatively estimate the accuracy of 
the tan /3 determination from measurements of the cross sections of these channels. For 
this purpose, we must first choose the relevant underlying parameters that fix the Higgs 
potential. In our analysis, we used a Higgs potential with only the leading correction 
from the top-stop loop, with stop mass 1 TeV. All the masses and mixings in the Higgs 
sector are then fixed if we determine the charged Higgs mass and tan {3. Notice that we 
use the simple form of the Higgs potential only to generate the event samples, and do 
not use any theoretical assumptions in the actual determination of tan j3. It is therefore 
straightforward to extend our analysis to the case with general radiative corrections. 

With the physical underlying parameters, we estimate the cross sections in channels 1 
- 8. These cross sections then determine the number of events that would actually 
be observed in each channel. We then postulate a hypothetical tan/?, calculate the 
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cross sections in channels 1-8 based on the postulated tan/?, and check whether the 
postulated tan/3 is consistent with observations. To quantify the argument, we define 



A*°= £ (S) 

CT -4- CT 

ixhannel stat 1 w syst 

where Ni (N[) is the number of events in channel i for the underlying (postulated) tan /?, 
and the quantities <7g tat and <7g yst are the statistical and systematic errors for channel i, 
respectively For simplicity, we add <7g tat and al yst in quadrature. The number of events 
has two origins: the signal and the background. Notice that the background cross section 
is independent of tan (3. Therefore, if some channel is dominated by the background, it 
cannot give a significant contribution to A% 2 . 

The statistical error is <7g tat = JN^ while the systematic error is given by 



•WW- (9) 

where the sum is over all quantities P that enter in the calculation of the numbers 
of events, and which therefore contribute systematic uncertainties. We include uncer- 
tainties from the following quantities (the uncertainties we use are in parentheses): the 
bottom quark mass (150 MeV |]J), cos 2 (/3 — a) (2 % [[!]]), the hadronic decay width of the 
heavy Higgses (20 % H), Br(H -> hh) (10 % 0), the 6-tagging efficiency (2 % H), 



and heavy Higgs masses (16 GeV /^0.035N H , where N H is the number of Higgs pair 
production events fT6|l). 

In the following, we will use Ax 2 to estimate the accuracy of the tan/3 measurement. 
We will show contours of constant A% 2 = 3.84, which we refer to as 95 % C.L. contours. 



4 NUMERICAL RESULTS 

In this section, we present a quantitative estimate of the expected uncertainty of the 
tan p measurement at the NLC Here, we assume the absence of SUSY decays modes, 
whose effects are discussed in the next section. 

We start with the NLC with = 500 GeV, and the charged Higgs mass = 200 
GeV. In this case, the decay of neutral Higgses into tt pair is kinematically forbidden. 
The tan f3 dependence of the cross sections then comes mainly from charged Higgs events. 
In Fig. Ill we plot the expected accuracy for tan/3 as a function of the input value (i.e., 
the actual value) of tan/3. Here, we use four different integrated luminosities: 25, 50, 
100, and 200 fb^ 1 , which correspond to 0.5, 1, 2, and 4 years of running with design 
luminosity. The qualitative behavior of the figure can be understood in the following 
way. If the actual value of tan/3 is close to 1, the branching ratio of the charged Higgs 
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decaying into tv t is too suppressed to be observed, and hence we may set an upper bound 
on tan/3 from the non-observation of the H ± leptonic decay mode. If the underlying 
tan/3 is in the moderate region (~ 3 - 10), we observe both the leptonic and hadronic 
decays of H ± . In this case, we can determine tan/3 from these modes, as can be seen in 
Fig. |l|. If the actual tan (3 is larger than ~ 10, pair production of heavy Higgses allows 
us to set a lower limit on tan/3. However, if tan/3 is large enough, the cross section for 
the process e + e~ — > tbH^ (channel 2) can be large enough to be observed, and both an 
upper and a lower limit on tan (3 can be obtained from channel 2. 

If the heavy Higgses are not kinematically accessible, or even if they are, it is advan- 
tageous to increase the beam energy. Therefore, we next present results for phase II of 
the NLC with ^ = 1 TeV. In Figs. g - 1, we present results for m H ± = 200 GeV, 300 
GeV, and 400 GeV, respectively. For the integrated luminosities, we use 100, 200, 400, 
and 800 fb -1 . In particular, for m#± = 200 GeV, we can see a great improvement of the 
result, comparing Figs. [I] andg. For m H ± = 300 GeV, we can still expect a precise deter- 
mination of tan /3, especially if the underlying tan j3 is moderate or large. For m#± = 400 
GeV, the result is noticeably worse. There are mainly two reasons for this. First, y/s is 
close to the threshold energy in this case. The Higgs production cross sections then be- 
comes small due to the phase space suppression, and the statistical errors become large. 
Furthermore, as the charged Higgs mass gets larger, the phase space suppression of the 
process H ± — > tb becomes less significant, and the decay mode H ± — > ru T is relatively 
suppressed. As a result, the branching ratio of the charged Higgs loses its sensitivity 
to tan/3, contrary to the cases with m#± = 200 and 300 GeV. In fact, in the case with 
rriH± = 400 GeV, tan/3 is mainly constrained by the AH production process with 46, 
2b2t, and 4t final states. Note, however, that even in this case, we may distinguish low, 
moderate, and high tan/3, and we can still obtain accurate measurements of tan/3 if the 
underlying value of tan j3 is in the moderate region. 

Before closing this section, we comment on the effects of the supersymmetric decay 
modes. Up to now, we have assumed that the Higgs bosons decay only into standard 
model particles. However, especially when the Higgs masses are large, Higgses may decay 
into superparticles, which one might think would lead to new sources of large systematic 
errors. Here, we would like to point out that this is not necessarily the case. Let us 
start with the decay mode H — > IrI* r . We can predict the branching ratio for this mode, 
since this process is induced by D-term interactions^. Thus, we only have to include this 
mode into the fit, and no new large systematic uncertainties arise. The primary effect 
of this decay mode is then to reduce the number of events from AH production, and in 
fact, the numerical results are not changed much ||. The same argument applies to the 
decay of H into pairs of left-chirality sfermion. However, if the heavy Higgses can decay 
into sfermion pairs with different chirality, such as H — > IrI* l , the branching ratios of the 

■^Slepton masses will be measured at the NLC with good accuracy. 
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heavy Higgses depend on additional MSSM parameters, such as fi and the trilinear A 
terms. These decay modes may then be a source of a large systematic errors. Of course, 
the new parameters may also be measured from different observables; for example, /x 
may be measured from chargino and neutralino masses, and the trilinear scalar couplings 
may be measured from left-right mixings. A complete analysis would therefore require 
a simultaneous fit to all of these parameters. 

Finally, we briefly consider decays to charginos and neutralinos. These decay may be 
dominant in some regions of parameter space. However, if only decays to the lighter two 
neutralinos and the lighter chargino are available, and these are either all gaugino-like 
or all Higgsino-like, as is often the case, these decays are suppressed by mixing angles. 
If we are in the mixed region, these decay rates may be large, but in this case, all six 
charginos and neutralinos should be produced, and the phenomenology is quite rich and 
complicated. 



5 DISCUSSION 



The results presented in the previous section indicate that tan (3 may be well determined 
from the study of production and decay of the heavy Higgses at the NLC The uncertainty 
in the measured tan/3 can be O(10 %) or less, depending on the underlying parameters. 
For example, for y/s = 500 GeV, m H ± = 200 GeV, and £ = 100 fb _1 , the observed 
value tan /3 b s will be in the following ranges for various input values of tan /3: 



tan (3 (input) = 2 
tan f3 (input) = 3 
tan (3 (input) = 5 
tan/? (input) = 10 
tan (3 (input) = 60 



tan/? obs < 2.9, 
2.5 < tan (3 ohs < 3.6, 

4.5 < tan /3 b s < 5.5, 

7.6 < tan j3 ohs < 30, 
40 < tan/3 obs < 90. 



In our study, we have not adopted any assumption which strongly depends on a 
specific model. Thus, if the heavy Higgses are kinematically accessible at the NLC, 
we can expect to obtain a constraint on tan (3 from the detailed study of heavy Higgs 
bosons. 

Finally, we discuss implications of precise determinations of tan/3. First of all, it 
should be emphasized that a determination of tan (3 may help us understand physics at 
very high energy scales. For example, the simple SO (10) GUT predicts large values of 
tan/3 |I7J. The unification of mf,-m T based on SU(5)-type GUTs is an another example. 
It prefers a value of tan/3 that is very large or close to 1 [18]. Precise determinations of 
tan (3 can be excellent tests of these scenarios. 

The parameter tan j3 is also important for the determination of the SUSY breaking 
scalar masses. Neglecting mixings, the physical masses of sfermions are the sum of soft 
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SUSY breaking masses and the D-term contribution which depends on tan (3. Thus, we 
must know tan (3 to determine soft SUSY breaking mass parameters from the physical 
masses of sfermions. If tan/? is completely unknown, soft SUSY breaking masses may 
have large uncertainties even if we can measure the sfermion masses very accurately. 
Measurement of tan/3 reduce this uncertainty. The important point is that the D- 
term contribution is proportional to cos 2/3, and hence it becomes insensitive to tan/3 
once tan/3 is larger than 3-4. Given our result, the uncertainty related to the D- 
term contribution is smaller than the error from the sfermion mass measurement, if 
tan (3 is larger than 3-4. The spectrum of the soft SUSY breaking masses may contain 
information about the origin of SUSY breaking [19] and the gauge and/or flavor structure 



at high scales [20 



Furthermore, if we combine the tan j3 determination with measurements of the Higgs 
masses, we may be able to gain some information about the Higgs potential. The 
Higgs potential is strongly constrained at the tree level, but radiative corrections are 
quite significant. In particular, the top squark plays an important role, and the precise 
determination of tan (3 as well as the Higgs masses may give us some constraint on top 
squark masses. 

In summary, determination of the tan f3 parameter can potentially give us rich infor- 
mation about physics at the high scale, including the possibility of GUTs, high energy 
flavor structures, and the origin of SUSY breaking. 
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Figure 1: Accuracy of the measured tan/5 (95 % C.L.) for y/s = 500 GeV, m H ± = 200 
GeV, eb = 60 %, and four integrated luminosities: 25, 50, 100, and 200 fb" 1 (from 
outside to inside). 
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Figure 2: Accuracy of the measured tan/5 (95 % C.L.) for ^Js = 1 TeV, m#± = 200 
GeV, e b = 60 %, and four integrated luminosities: 100, 200, 400, and 800 fb" 1 (from 
outside to inside). 
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Figure 3: Same as Fig. 0, except for m H ± = 300 GeV. 
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Figure 4: Same as Fig. ^ except for m#± = 400 GeV. 
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